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Abstract  -  Track-Before-Detect  (TkBD)  is  a  paradigm 
that  combines  the  target  detection  and  estimation  processes 
that  are  usually  sequentially  applied  to  sensor  data  in  a  con¬ 
ventional  system.  Under  TkBD  the  single  frame  detector  is 
removed  and  the  tracker  is  supplied  with  the  whole  sensor 
image.  Detection  decisions  are  then  shifted  to  the  output 
of  the  tracker  which  is  able  to  use  temporal  correlation  to 
improve  the  decision  performance.  A  fundamental  measure 
used  by  most  TkBD  approaches  is  the  likelihood  ratio  of  the 
sensor  data  and  this  is  formed  as  the  product  over  individual 
cell  likelihoods  under  the  assumption  of  spatially  indepen¬ 
dent  noise.  However,  that  approach  exploits  only  the  enve¬ 
lope  of  the  known  sensor  point  spread  function.  This  article 
presents  an  approach  for  determining  the  data  likelihood  ra¬ 
tio  that  also  includes  phase  information.  This  alternative 
likelihood  ratio  formulation  is  shown  to  both  improve  the 
discrimination  of  targets  from  noise  and  reduce  the  compu¬ 
tation  overhead  of  the  algorithm. 

1  Introduction 

Traditional  tracking  algorithms  are  designed  assuming 
that  the  sensor  provides  a  set  of  point  measurements  at  each 
scan.  The  tracking  algorithm  links  measurements  across 
time  and  estimates  parameters  of  interest.  However,  a  practi¬ 
cal  sensor  may  provide  a  data  image,  where  each  pixel  corre¬ 
sponds  to  the  received  power  in  a  particular  spatial  location 
(e.g.  range  bins  and  azimuth  beams).  In  this  case,  the  com¬ 
mon  approach  is  to  apply  a  threshold  to  the  data  and  to  treat 
those  cells  that  exceed  the  threshold  as  point  measurements, 
perhaps  using  interpolation  methods  to  improve  accuracy. 
This  is  acceptable  if  the  signal  to  noise  ratio  (SNR)  is  high. 
For  low  SNR  targets  the  threshold  must  be  low  to  allow  suf¬ 
ficient  probability  of  target  detection.  A  low  threshold  also 
gives  a  high  rate  of  false  detections,  which  causes  the  tracker 
to  form  false  tracks.  An  alternative  approach,  referred  to  as 
track-before-detect  (TkBD),  is  to  supply  the  tracker  with  all 
of  the  sensor  data  without  applying  a  threshold.  This  im¬ 
proves  track  accuracy  and  allows  the  tracker  to  follow  low 
SNR  targets. 
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The  key  difference  between  TkBD  algorithms  and  con¬ 
ventional  trackers  is  the  measurement  model.  Conventional 
trackers  use  a  point  measurement  model  that  typically  as¬ 
sumes  that  the  observations  are  kinematic  states  of  the  tar¬ 
get  corrupted  with  additive  (usually  Gaussian)  noise  (e.g. 
[1,  3]).  However,  a  TkBD  algorithm  requires  a  statistical 
model  of  the  signal  received  in  each  sensor  pixel. 

Earlier  TkBD  efforts  focussed  on  optical  sensors  [2,  6,  7, 
13]  and  assumed  that  the  sensor  noise  was  Gaussian.  More 
recent  efforts  have  considered  the  radar  problem  where  the 
pixel  intensity  is  either  Ricean  or  Rayleigh  distributed  de¬ 
pending  on  the  presence  or  absence  of  targets  [4,  9,  10,  11]. 
The  overall  likelihood  ratio  is  formed  by  taking  a  product 
of  the  pixel  likelihood  ratios  for  all  pixels.  This  product  is 
a  result  of  assuming  that  the  noise  is  spatially  uncorrelated. 
The  Ricean  and  Rayleigh  distributions  are  functions  of  only 
the  envelope  of  the  sensor  data,  so  this  approach  does  not 
use  phase  information. 

This  article  demonstrates  an  alternative  likelihood  formu¬ 
lation  that  allows  for  spatially  correlated  noise  and  includes 
the  phase  information.  The  performance  of  a  grid  based  re¬ 
cursive  posterior  filter  [5,  14]  is  analysed  for  the  traditional 
likelihood  function  using  only  envelope  data  and  the  alter¬ 
native  likelihood  function  that  includes  phase.  The  inclu¬ 
sion  of  phase  information  is  shown  to  improve  the  detection 
and  estimation  performance  of  the  filter.  Furthermore,  the 
complex  likelihood  function  is  an  order  of  magnitude  more 
efficient  to  evaluate  than  the  envelope  one. 


2  Problem  Definition 

As  in  [9,  Ch.  11],  consider  a  sensor  that  collects  a  se¬ 
quence  of  two-dimensional  images  (frames).  When  present, 
a  target  moves  in  the  plane  according  to  a  known  statistical 
process.  Each  pixel  in  a  frame  is  a  complex  valued  random 
variable,  and  the  target  contribution  is  via  a  known  complex 
point  spread  function  with  an  arbitrary  and  unknown  phase 
offset. 
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2.1  Target  Model 

For  simplicity  of  notation,  assume  a  discrete  time  model, 
with  a  fixed  frame  period,  T.  The  target  state  at  frame  k, 
Xfc,  consists  of  position  and  velocity  in  the  plane,  and  the 
intensity  of  the  returned  signal,  Ik,  i.e. 

Xfc  =  [xk  Xk  yk  yk  4]T-  (1) 


The  evolution  of  the  state  is  modelled  by  the  linear  stochas¬ 
tic  process 

Xfc  =  Fxfc_i  +  Vfc,  (2) 


where  vk  is  a  noise  process  and  the  transition  matrix  is  given 
by 


F  = 
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0 
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(3) 


The  process  noise,  vk,  is  the  usual  Gaussian  random  vari¬ 
able  with  covariance  Q  given  by 
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where  qs  is  the  power  spectral  density  of  the  acceleration 
noise  in  the  spatial  dimensions  and  q,  is  the  power  spectral 
density  of  the  noise  in  the  rate  of  change  of  target  return 
intensity. 


2.3  Envelope  Likelihood 

The  standard  treatment  of  this  measurement  model  is  to 
assume  that  the  noise  is  spatially  uncorrelated  and  to  form 
the  likelihood  of  its  envelope  (i.e,  the  magnitude  of  the  data) 
[11,  4,  9],  Thus  R  is  assumed  to  be  an  identity  matrix  scaled 
by  a2. 

The  envelope  likelihood  considers  only  the  magnitude  of 
the  sensor  response  at  each  pixel.  The  pixel  magnitude  will 
be  Ricean  distributed,  if  there  is  a  target  present,  or  Rayleigh 
distributed  if  there  is  no  target  [12].  Thus  the  pdf  of  \zk  |  is 


V  \zk\ 


Xfc  = 


exp 


141  +|^(xfc)| 

2a2 

l^(xfc)  41 


X  In 


,  (6) 


if  the  target  is  present,  or 


P\  zk  = 


exp 


2a2 


(7) 


if  there  is  no  target,  where  Iq(-)  is  the  modified  Bessel  func¬ 
tion. 

Note  that  the  measurement  probability  implied  by  (5)  de¬ 
pends  on  the  phase  shift,  </>.  The  pdf  (6)  is  the  marginal 
distribution  achieved  by  the  integral 


2.2  Measurement  Model 

The  measurement  at  each  time  is  an  image  of  an  arbitrary 
dimensionality  with  N  pixels  indexed  by  i  =  1 ...  AT.  De¬ 
note  the  complex  response  in  the  ith  pixel  of  the  fcth  frame  as 
zlk  and  let  zk  be  a  stacked  vector  of  all  the  pixel  responses 
of  the  fcth  frame.  The  target  signal  at  the  ith  pixel  is  de¬ 
noted  by  /i'(xfc)  and  h(xfc)  is  a  stacked  vector  of  the  target 
signal  over  all  of  the  pixels.  This  target  signal  may  be  the 
point  spread  function  (psf)  of  the  sensor,  the  target  signa¬ 
ture  (range  profile)  or  a  combination  of  both.  It  is  shifted 
to  be  centred  over  the  target  state  and  scaled  by  the  target 
amplitude  and  assumed  to  be  known. 

Assuming  additive  noise,  the  measurement  equation  is 


P  (|4 

The  Bessel  function  arises  as  a  result  of  this  integral  [8], 
The  likelihood  ratio  for  the  magnitude  response  at  pixel  i 
is  thus 


Xfc  = 


Xfc,</>)f>(0)#.  (8) 


(9) 


zk  =  exp{j^}h(xfc)  +  rifc,  (5) 

where  <f>  is  an  unknown  phase  shift,  uniformly  distributed 
over  [0,  27t],  and  rifc  is  a  stacked  vector  of  the  (complex) 
pixel  noise  signals.  Assume  that  rifc  is  complex  Gaussian 
with  a  known  covariance  matrix,  R.  For  this  paper,  the  noise 
will  be  assumed  to  be  spatially  independent,  so  R  is  an  iden¬ 
tity  matrix  scaled  by  the  known  noise  variance,  a2 . 

The  target  peak  SNR  quantifies  the  height  of  the  peak  of 
the  target  point  spread  function  relative  to  the  noise  floor, 
and  represents  a  measure  of  how  easy  it  is  to  detect  the  tar¬ 
get.  The  point  spread  function  is  assumed  to  be  normalised 
such  that  the  contribution  to  cell  i  is  Ik  when  the  target  is 
located  exactly  on  the  sample  point  for  the  cell.  Thus  the 
peak  SNR  in  dB  is  given  by  20  log  { Ik  /cr2  } . 


Since  the  pixels  are  assumed  to  be  conditionally  indepen¬ 
dent,  the  likelihood  of  the  envelope  of  the  whole  image  is 
simply  the  product  over  the  pixels 


where  H  denotes  the  Hermition  conjugate. 

Typically  this  product  is  truncated  to  only  include  those 
pixels  that  have  a  significant  target  contribution,  hl  (xk),  to 
avoid  unnecessary  computation  [11]. 


619 


If  a  prior  distribution  is  assumed  for  the  target  intensity, 
p{Ik),  then  an  intensity  independent  marginal  likelihood  is 
given  by 


Ce{zk  |xfc)  =  /  Ce{zk\xk)p(Ik)dIk.  (11) 
Jo 

This  integral  may  be  numerically  approximated  using  a  sum¬ 
mation. 

An  important  implication  of  the  above  method  is  that 
marginalisation  over  the  random  phase  shift,  <f>,  has  been  re¬ 
peated  independently  for  each  pixel.  This  is  equivalent  to 
assuming  that  there  is  a  different  phase  shift  for  every  pixel, 
(f>\.  This  is  not  the  case,  and  so  the  envelope  likelihood  has 
not  captured  all  of  the  available  information. 

3  Complex  Likelihood 

The  factorised  form  of  the  envelope  likelihood  results  in 
a  loss  of  information  because  it  fails  to  include  the  phase  of 
the  sensor  data.  It  also  requires  the  calculation  of  a  Bessel 
function  for  each  pixel  with  a  significant  psf  contribution. 
Experience  shows  that  the  evaluation  of  these  Bessel  func¬ 
tions  is  by  far  the  most  computationally  expensive  part  of 
TkBD  algorithms  that  use  this  likelihood. 

Consider  now  the  joint  likelihood  of  the  whole  sensor  im¬ 
age.  For  compactness,  the  notation  will  omit  the  dependence 
of  the  psf  on  the  target  state,  omit  the  frame  index  k,  and  let 
s  =  exp{j</>}  with  s *  its  conjugate.  This  likelihood  defines 
a  matched  filter  which  can  account  for  both  spatially  corre¬ 
lated  noise  and  the  phase  response  of  the  psf.  The  derivation 
of  this  likelihood  is  now  presented. 

The  probability  density  of  z  when  a  target  is  present  is 

P(z|target,  <j>)  = 

— — rexp(-i  (z  _  sh)H  FT1  (z  -  sh)}  ,  (12) 

|2ttR|3  l  2  J 

and  when  a  target  is  absent  is 

P(z|no  target)  =  - r  exp  j  — -zHR_1z  1  ,  (13) 

|2ttR|3  l  2  J 

since  n  is  Gaussian  with  covariance  R. 

Thus  the  likelihood  ratio  is 

C{z\x,4>) 

=  exp  |  —  }  (z  —  S‘h)H  R-1  (z  —  sh)  +  }zHR_1zl , 


where 


A  =  exp|-^hHR  ^j, 

£  =  hHR_1z  =  2exp{j<9}. 


Algebraic  manipulation  yields 
C(z\x,4>) 

=  A  exp  |  }  (cos  <j)  +  j  sin  </>)£*  +  }  (cos  <j>  —  j  sin  <f>)£  , 
I  z  z 

=  A  exP  |  \  cos  (j)[C  +  £]  +  p  sin  </>[£*  -  £] 


=  A  exp  <  2  cos  (j>  cos  9  +  S  sin  <f)  sin  ( 


=  A  exp  <  Ecos(<j>  —  9)  >. 


(15) 


Marginalise  over  (j>  assuming  a  uniform  distribution: 


C{z\x) 


C(z\x ,  (f))p{(j))d(j) 


A  rn  (  ) 

—  J  exp  |scos(</>  —  9)  jd(f> 


exp  \  — -hHR-1h  \  l0  (  |hHR_1z| 


If  the  noise  is  spatially  uncorrelated  with  variance  a2, 
then  R  is  the  identity  matrix  scaled  by  a2  and  the  likelihood 
simplifies  to 


C{z\x) 


(17) 


As  with  the  envelope  likelihood,  the  above  expression  is 
implicitly  dependent  on  the  target  amplitude.  One  way  to 
remove  this  dependence  is  to  marginalise  over  it  as  in  (11). 

The  likelihood  expression  is  very  similar  to  the  envelope 
likelihood  except  that  whereas  the  envelope  likelihood  in¬ 
cludes  a  product  of  Bessel  functions,  (17)  has  a  single  Bessel 
function  evaluated  at  the  output  of  a  matched  filter.  Because 
evaluating  the  Bessel  functions  is  the  main  computation  cost 
of  many  TkBD  algorithms,  one  would  expect  a  significant 
speed  improvement  by  using  (17)  instead  of  the  envelope 
likelihood.  The  key  difference  is  that  the  matched  filter  form 
includes  the  phase  profile  of  the  psf  which  is  lost  in  the  en¬ 
velope  likelihood.  This  extra  information  should  lead  to  su¬ 
perior  performance. 


=  exp  | ^zHR-1sh  +  is*hHR-1z  -  ^s*hHR"1sh| , 

s*hHR“1z|  , 

=  Aexpj^C  +  <14) 


=  AexpjiszHR  1h+i 


4  State  Posterior  Recursion 

The  test  TkBD  algorithm  is  a  discretised  numerical  ap¬ 
proximation  to  the  exact  Bayesian  posterior  recursion.  The 
algorithm  uses  a  fixed  grid  in  the  state  space  and  propagates 
the  posterior  probability  mass  of  the  target  state  over  this 
grid.  The  mode  of  the  posterior  is  used  as  the  estimate  in 
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this  paper,  since  it  is  likely  that  the  distribution  will  be  multi¬ 
modal,  and  the  mean  may  be  sensitive  to  distant  low  likeli¬ 
hood  modes. 

The  posterior  pdf  of  the  target  state  can  be  recursively 
determined  using  the  well  known  Bayesian  relationship 


p(xk\Zk)  oc  p(zk\xk)  /  p(xfc|xfc_i)p(xfc_i|Zfe_i)dxfc_i, 

(18) 

where  Zk  =  z\...  zk- 

The  estimator  in  this  article  uses  a  direct  approximation 
to  (18)  based  on  a  discretisation  of  the  state  space.  Choose 
a  uniformly  spaced  set  of  states,  X  (which  is  not  necessarily 
related  to  the  discrete  measurement  function)  such  that 


xk(q,r,s,t)  = 


A  xq  Ays  =ff  Ik 


(19) 


for  some  integers  q,r,s  and  t. 

(18)  can  then  be  approximated  by 

p(xk\Zk)  « KjC(zk\xk)  ^  p(xfc|xfc_i)p(xfc_i|Zfc_i), 

(20) 

where  K  is  a  normalising  constant  and 


F[x{q,r,s,t)). 

xGrV  q=—oor=—oos=—oot=—oo  '  ' 

(21) 

The  approximation  is  exact  in  the  limit  as  X  approaches 
•ft4  for  non-negative  continuous  functions,  as  used  here. 

The  first  term  in  (20)  is  the  intensity  independent  marginal 
likelihood,  defined  by  (11)  and  using  either  the  envelope 
likelihood  or  the  matched  likelihood. 

The  state  evolution  term  in  (20),  p(xfc|xfc_i),  is  a  four 
dimensional  matrix.  The  dynamics  determines  the  centre  of 
the  distribution,  i.e.  Fxfc_ls  and  the  shape  is  prescribed  by 
the  probability  mass  of  the  process  noise,  vk,  which  is  now 
discrete  on  X.  In  order  to  minimise  the  computation  burden 
of  (20),  the  probability  mass  of  v/i;  is  assumed  to  be  zero 
outside  a  compact  region  centred  on  the  origin. 

The  discrete  state  space  is  augmented  with  a  null  state, 
0,  to  indicate  the  possibility  that  there  is  no  target.  Denote 
the  probability  of  target  death  as  P,i,  and  the  probability  of 
target  birth  as  Then  the  evolution  probability  in  (20)  is 
given  by 


p(xfc|xfc_i) 


{1  -  Pb  xfc  =  0,xfe_ i=0, 

Pd  xfc  =  0,xfe_  i^0, 

Pb/\X\  Xfc  ^  0,Xfc_i  =  0, 

(1  -Pd)Pv{Vk)  xfc^0,xfc_!^0, 

(22) 


where  \X\  is  the  number  of  discrete  states  in  X,  vk  = 
xk  —  Fxk-i  and  pv(vk)  is  the  process  noise  probability  dis¬ 
tribution  (refer  to  section  2.1). 

The  parameters  Pb  and  Pd  are  tuning  parameters  and  may 
be  adjusted  to  optimise  detection  performance.  The  selec¬ 
tion  of  the  state  space,  X,  is  a  trade-off  between  estima¬ 
tion  accuracy,  which  improves  with  finer  resolution,  and 


computation  requirement,  which  increases  with  \X\.  The 
process  noise  pmf  also  affects  estimation  accuracy,  as  well 
as  providing  some  capacity  to  handle  model  mismatch  be¬ 
tween  the  assumed  target  model  and  the  true  target  mo¬ 
tion.  The  algorithm  is  initialised  with  p(x o  =  0)  =  1  and 
p(x0)  =  0  V  x0  /  0. 

Once  the  pdf  of  the  state  has  been  evaluated,  a  state  esti¬ 
mate  can  be  obtained  by  selecting  the  state  with  the  highest 
probability.  In  the  event  that  this  state  is  the  null  state,  then 
the  algorithm  reports  that  there  is  no  target.  To  account  for 
the  case  where  the  pdf  has  a  peak  that  is  spread  over  several 
grid  cells,  the  implementation  used  in  this  article  finds  the 
highest  probability  non-null  state  and  accumulates  the  prob¬ 
ability  in  the  adjacent  cells.  If  the  accumulated  probability 
is  higher  than  the  null  state  probability,  then  a  detection  is 
reported. 

In  the  experimental  analysis,  it  will  be  seen  that  imple¬ 
mentations  of  the  Bayesian  posterior  recursion  using  the  two 
likelihood  alternatives  have  substantially  different  computa¬ 
tion  costs.  In  order  to  remove  the  effects  of  implementation 
efficiency,  the  computation  cost  of  the  algorithm  is  now  con¬ 
sidered  analytically. 

Assume  that  the  state  space  X  is  a  regular  grid  such  that 
\X\  =  Ax  Ay,  where  Ax  is  the  number  of  position  states 
in  X  and  Ay  is  the  number  of  velocity  states.  Let  V  denote 
the  number  of  elements  of  X  for  which  the  process  noise  has 
a  non-zero  probability  mass. 

The  likelihood  is  independent  of  the  target  speed  and  so 
an  efficient  implementation  requires  Ax  likelihood  calcula¬ 
tions,  each  of  which  has  a  cost  Cc- 

For  each  state  locale,  (20)  requires  V  products  and  accu¬ 
mulates  under  the  sum  with  a  further  two  multiplies  outside 
the  sum.  Determining  the  scaling  constant  requires  a  sum 
over  the  whole  state  space.  Let  Cx  and  CL  denote  the  cost 
of  a  multiplication  and  a  sum  respectively. 

Thus  the  total  cost  of  a  single  time  slice  of  the  estimator 
is 


C  =  AxCc  +  AxAaC+ 

+  ^x  Ay  +  2)CX  +  V  C'+ 1 , 

~  AxC c  +  AxAjrV  (Cx  +  C+) ,  (23) 

since  V  1. 

If  the  state  space  has  a  fixed  physical  extent  in  space  and 
velocity  and  a  symmetrical  sampling  spatial  frequency  fs 
then  the  cost  can  be  written  as 

C  »  tfsfaiCc  +  (fe)8a2  (Cx  +  C+) ,  (24) 

where  a±  and  a 2  are  constants  whose  value  are  a  function  of 
the  physical  extent  of  the  state  space  in  position  and  velocity. 

So  the  algorithm  cost  contains  a  term  that  has  0(N2) 
complexity  in  the  grid  resolution  due  to  the  likelihoods  and  a 
term  that  has  0(N8)  complexity  due  to  the  dynamic  model. 

It  is  intuitive  that  the  cost  of  calculating  the  envelope 
likelihood  is  much  higher  than  that  of  the  matched  likeli¬ 
hood.  Although  the  second  term  is  high  in  complexity,  it  is 
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amenable  to  pipeline  processing  and  may  be  relatively  cheap 
for  problems  over  a  small  spatial  area. 

5  Performance  Evaluation 

Bayesian  algorithms  constructed  using  the  two  likelihood 
functions  are  now  compared  via  Monte  Carlo  simulation. 
The  test  scenarios  consist  of  a  small  surveillance  volume 
containing  either  a  single  straight  line  target,  or  no  target. 
Although  the  straight  line  target  assumption  is  rather  limit¬ 
ing,  this  paper  is  focussed  on  the  measurement  model,  not 
the  target  dynamics,  so  there  is  no  requirement  for  more 
challenging  targets. 

A  TkBD  algorithm  is  primarily  a  detector  that  exploits 
a  target  dynamics  model  to  achieve  integration  gain.  There¬ 
fore  the  most  important  measure  of  performance  is  how  well 
the  algorithm  discriminates  between  the  presence  and  the 
absence  of  a  target.  Of  secondary  interest  is  the  processing 
resource  required  to  employ  the  algorithm.  This  quantity  is 
highly  dependent  on  the  implementation.  However,  the  vari¬ 
ation  in  implementation  cost  should  reflect  the  complexity 
model  described  above.  The  estimation  error  is  also  of  in¬ 
terest,  although  it  is  much  less  important  than  the  detection 
performance  and  processing  cost. 

5.1  Simulation  Details 

The  simulation  scenario  uses  a  30  x  45  pixel  scan  with  unit 
variance  complex  normal  noise.  When  present,  the  target 
signature  is  the  spectral  response  of  a  Hann  window  centred 
at  the  target  location.  The  target  moves  with  a  fixed  speed  of 
1  pixels  per  frame.  For  each  Monte  Carlo  trial  a  heading  is 
uniformly  sampled  between  0  degrees  (East)  and  45  degrees 
(North-East). 

The  target  state  space  was  defined  as  the  measurement 
domain  in  the  spatial  dimensions  and  ±3  pixels  per  frame  in 
each  velocity  dimension.  The  resulting  4D  state  space  has  a 
size  of  30  x  7  x  45  x  7. 

Scalloping  loss  may  reduce  the  sensitivity  of  the  algo¬ 
rithm  to  a  target  that  is  located  part  way  between  sample 
points,  resulting  in  a  potential  variation  in  performance  de¬ 
pending  on  the  relative  position  on  the  target  with  respect 
to  the  grid.  In  order  to  average  over  this  potential  variation, 
the  initial  target  position  for  each  Monte  Carlo  trial  was  ran¬ 
domly  sampled  from  the  range  2.5  ...  3  independently  in  X 
and  Y. 

The  Bayesian  system  model  comprises  a  target  dynam¬ 
ics  model,  a  measurement  model,  and  a  target  birth-death 
model.  Of  these,  the  measurement  model  is  assumed  to  be 
known,  and  the  target  model  is  assumed  to  be  fixed.  The 
target  model  is  specified  by  the  process  noise,  which  is  as¬ 
sumed  zero  outside  ofa3x3x3x3  region  with  probability 
0.7  at  the  origin  and  the  remaining  0.3  spread  equally  over 
the  other  80  elements. 

The  parameters  of  the  birth-death  model  are  the  probabil¬ 
ity  of  new  target  birth  and  the  probability  of  old  target  death. 
For  this  comparison,  it  was  found  that  using  a  very  small 
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Figure  1 :  Overall  ROC 

probability  of  death  gave  best  performance  for  all  birth  prob¬ 
abilities,  so  the  probability  of  death  was  fixed  at  Pd  =  10-5. 

The  probability  of  birth,  /  ’,,  was  varied  to  trace  out  a  Re¬ 
ceiver  Operating  Characteristic  (ROC)  curve.  For  each  value 
of  Pb,  100  Monte  Carlo  trials  were  performed  with  a  target 
present  and  100  without  a  target. 

5.2  Performance  Measures  and  Results 

Five  Measures  Of  Performance  (MOPs)  were  defined  for 
the  comparison:  overall  ROC;  per-scan  ROC;  per-scan  de¬ 
tection  as  a  function  of  SNR;  RMS  position  error  and  com¬ 
putation  resource. 

For  TkBD,  the  main  focus  is  on  detection  of  targets,  thus 
the  first  three  measures  are  the  most  important.  Of  lesser  in¬ 
terest  is  the  estimation  accuracy  and  computation  cost,  pro¬ 
vided  that  the  latter  is  not  unreasonably  high. 

When  a  target  scenario  was  tested,  the  target  was  declared 
detected  on  a  particular  scan  if  the  TkBD  algorithm  pro¬ 
duced  an  output  track  for  that  scan  that  was  sufficiently  close 
to  the  true  target  location.  The  association  gate  used  only  the 
target  position  and  had  a  radius  of  2  pixels. 

5.2.1  Overall  ROC 

The  proportion  of  targets  detected  is  defined  as  the  frac¬ 
tion  of  Monte  Carlo  trials  where  the  TkBD  algorithm  pro¬ 
duced  a  track  that  was  within  the  association  gate  of  the  tar¬ 
get  location  for  at  least  one  scan.  The  proportion  of  trials 
with  false  tracks  is  defined  as  the  fraction  of  Monte  Carlo 
trials  where  the  TkBD  algorithm  produced  at  least  one  false 
track.  This  statistic  only  includes  trials  where  there  was  no 
target  present. 

The  overall  ROC  was  generated  by  plotting  the  proportion 
of  targets  detected  as  a  function  of  the  proportion  of  trials 
with  false  tracks.  This  MOP  quantifies  how  well  the  differ¬ 
ent  likelihood  ratios  discriminate  between  data  that  contains 
a  target  and  data  that  does  not. 

Figure  1  shows  the  overall  ROC  curves  for  both  of  the 
likelihood  functions  at  target  SNR  values  of  6dB  and  3dB. 
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Figure  2:  Per-scan  ROC 

Both  likelihood  functions  detect  all  targets  at  6dB,  even  for 
very  low  false  alarm  rates.  At  3dB,  the  matched  likelihood 
still  provides  very  good  detection  performance,  detecting 
approximately  95%  of  targets  at  very  low  false  alarm.  How¬ 
ever,  the  envelope  likelihood  is  unreliable  at  3dB  and  only 
detects  around  half  of  the  targets.  The  matched  likelihood 
clearly  provides  superior  overall  detection  performance. 

5.2.2  Per-scan  ROC 

The  overall  ROC  does  not  include  any  information  about 
how  long  the  algorithm  takes  to  detect  a  target,  or  how  long 
false  tracks  persist.  In  order  to  limit  the  number  of  MOPs 
these  considerations  are  summarised  by  empirically  estimat¬ 
ing  the  per-scan  probability  of  detection  and  false  alarm. 

The  per-scan  proportion  of  targets  detected  is  defined  as 
the  fraction  of  scans  where  the  TkBD  algorithm  produced  a 
track  that  was  within  the  association  gate  of  the  target  loca¬ 
tion.  An  algorithm  that  quickly  initiates  a  track  and  main¬ 
tains  the  track  will  have  a  high  per-scan  proportion  detected, 
whereas  one  that  has  a  long  initiation  delay  or  quickly  di¬ 
verges  will  have  a  low  one. 

The  per-scan  false  track  rate  is  defined  as  the  average 
number  of  false  tracks  reported  by  the  TkBD  algorithm  per 
scan.  This  statistic  summarises  both  the  number  and  length 
of  false  tracks  and  is  simply  the  ratio  of  the  total  number  of 
false  reports  from  the  TkBD  output  to  the  total  number  of 
simulated  frames. 

The  per-scan  ROC  was  generated  by  plotting  the  per-scan 
proportion  detected  as  a  function  of  the  per-scan  false  track 
rate  and  is  shown  in  figure  2. 

At  6dB,  the  matched  likelihood  per-scan  proportion  de¬ 
tected  is  close  to  0.95,  reflecting  that  the  algorithm  estab¬ 
lishes  a  track  (on  average)  by  the  second  scan.  The  envelope 
likelihood  has  a  slightly  lower  per-scan  proportion  detected, 
reflecting  a  track  establishment  delay  of  around  2-3  scans  at 
6dB. 

At  3dB,  the  matched  likelihood  per-scan  proportion  de¬ 
tected  is  significantly  lower  than  the  overall  proportion  de- 
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Figure  3:  Per-scan  detection  as  a  function  of  SNR 


tected.  This  indicates  that  the  algorithm  now  takes  a  more 
substantial  time  to  establish  track,  on  average  8  scans  at  low 
false  alarm  rates.  The  envelope  likelihood  per-scan  propor¬ 
tion  detected  is  less  than  half  the  overall  proportion  detected. 
This  is  because  many  of  the  targets  it  detects  produce  only 
short  track  segments. 

As  before,  the  matched  likelihood  is  clearly  superior. 

5.2.3  Per-scan  detection  as  a  function  of  SNR 

The  ROC  curves  shown  in  figures  1  and  2  correspond  to 
two  SNR  values.  For  this  MOP,  the  per-scan  false  track  rate 
was  held  fixed  and  the  per-scan  proportion  of  targets  de¬ 
tected  was  observed  as  a  function  of  SNR.  For  this  MOP 
the  birth  probability  was  fixed  to  achieve  0.05  false  reports 
per  scan. 

The  per-scan  proportion  detected  is  shown  in  figure  3. 
It  demonstrates  that  the  improved  detection  ability  of  the 
matched  likelihood  function  is  approximately  equivalent  to 
a  gain  of  1.5  dB  in  signal  power. 

5.2.4  RMS  position  error 

The  two  dimensional  position  error  was  averaged  over 
those  frames  when  the  target  was  detected.  Figure  4  shows 
the  RMS  error  for  each  of  the  likelihoods  as  a  function  of 
SNR.  The  probability  of  birth  was  again  fixed  to  achieve 
0.05  false  reports  per  scan.  The  RMS  error  for  the  envelope 
likelihood  was  only  measured  at  3  dB  and  above  since  there 
are  too  few  tracks  to  get  reliable  estimates  below  this. 

As  should  be  expected,  the  estimation  error  is  lower  when 
the  target  SNR  is  higher.  The  RMS  error  of  the  matched 
likelihood  tracks  is  roughly  equivalent  to  that  of  the  enve¬ 
lope  likelihood  at  2  dB  higher  SNR. 

5.2.5  Computation  resource 

Table  1  gives  the  computation  resource  spent  by  the  two 
alternatives  both  in  absolute  CPU  seconds  and  in  relative 
terms.  Profiling  the  algorithm  code  using  the  envelope  like¬ 
lihood  function  showed  that  more  than  90  percent  of  the  cpu 
time  was  spent  in  the  likelihood  calculation,  at  around  55 
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Figure  4:  RMS  error  (pixels)  as  a  function  of  SNR 


likelihood 

cpu  time  (s) 

cpu  time  (ratio) 

envelope 

39,100 

11 

matched 

3,570 

1 

Table  1 :  Computation  requirements 


seconds  per  monte  carlo  trial.  In  contrast,  the  calculations 
for  the  matched  likelihood  took  only  1.4  seconds  per  trial, 
on  average.  This  reduction  is  partly  due  to  the  number  of 
Bessel  function  evaluations,  which  dropped  by  a  factor  of 
approximately  20.  However,  additional  gain  was  made  due 
to  simplification  of  the  inner  loop  in  the  Matlab  code.  This 
additional  gain  is  not  likely  to  be  reflected  in  an  implementa¬ 
tion  of  the  algorithm  using  a  lower  level  software  language. 
The  nett  result  is  approximately  an  order  in  magnitude  speed 
improvement  in  the  total  algorithm. 

The  difference  in  computation  resource  required  for  an 
operational  implementation  of  these  algorithms  would  be 
highly  dependent  on  how  the  Bessel  functions  were  evalu¬ 
ated.  If  low-accuracy  approximations  to  the  likelihood  were 
found  to  be  sufficient,  then  the  speed  difference  may  be 
small.  However,  if  higher  accuracy  Bessel  evaluations  were 
required,  the  order  of  magnitude  speed  difference  may  be  re¬ 
alistic,  since  Matlab  uses  external  optimised  code  to  perform 
this  task. 

6  Summary 

This  paper  demonstrated  that  the  performance  of  a 
Bayesian  Track-Before-Detect  algorithm  was  significantly 
improved  by  making  use  of  the  phase  signature  of  the  target 
in  addition  to  the  amplitude  signature.  Intuitively,  this  phase 
information  may  be  incorporated  by  using  a  matched  filter 
tuned  to  the  known  complete  point  spread  function  of  the 
sensor.  This  matched  filter  is  in  addition  to  those  that  may 
be  used  for  range  gating  or  other  pre-tracker  processing. 


The  use  of  phase  information  was  found  to  be  roughly 
equivalent  to  a  1.5  dB  improvement  in  signal  power,  com¬ 
pared  with  the  phase  ignorant  algorithm.  The  phase  infor¬ 
mation  not  only  improved  detection  and  estimation  perfor¬ 
mance,  but  also  lead  to  a  more  efficient  algorithm  with  an  or¬ 
der  of  magnitude  lower  computation  resource  requirements. 
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